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We present measurements of the branching fraction and time-dependent CP asymmetries in B° — ► 
J/V>7T° decays based on 466 million T(4S) — > BB events collected with the BABAR detector at the 
SLAC PEP-II asymmetric-energy B factory. We measure the CP asymmetry parameters S = 
-1.23 ± 0.21(stat) ± 0.04(syst) and C = -0.20 ± 0.19(stat) ± 0.03(syst), where the measured value 
of S is 4.0 standard deviations from zero including systematic uncertainties. The branching fraction 
is determined to be B(B° -> J/ipn ) = (1.69 ± 0.14(stat) ± 0.07(syst)) x 10~ 5 . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charge conjugation-parity (CP) violation in the B meson system has been established by the BABAR [1] 



4 



and Belle [2] collaborations. The Standard Model (SM) 
of electroweak interactions describes CP violation as a 
consequence of a complex phase in the three-generation 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing ma- 
trix [3] . Measurements of CP asymmetries in the proper- 
time distribution of neutral B decays to CP eigenstates 
containing a J/ip and K° meson provide a precise mea- 
surement of sin 2(3 [4], where (3 is arg [ -V cd V* h /V td V t ^ } 
and the are CKM matrix elements with i, j quark 
indices. 

The decay B° — > J/ipir° is a Cabibbo-suppressed 
b — > ccd transition to a CP-even final state whose tree 
amplitude has the same weak phase as the b — -> ccs 
modes, e.g., the decay B° — > J/ipK®. The b — > ccd 
loop (penguin) amplitudes have different weak phases 
than the tree amplitude. If there is a significant pen- 
guin amplitude in B — > J/tpir°, then the measured val- 
ues of the CP asymmetry coefficients S and C will differ 
from the tree level expectations of —sin 2/3 and 0, re- 
spectively, and this mode could be sensitive to physics 
beyond the SM [5]. The coefficient 5" is related to CP 
violation in interference between amplitudes of direct de- 
cay, and decay after mixing, and C is related to direct 
CP violation. An additional motivation for measuring 
S and C from B° — > J/ipw is that they can provide a 
model-independent constraint on the penguin contami- 
nation within B a -» J/ipK° s [6]. 

The data used in this analysis were collected with the 
BABAR detector [7] at the PEP-II asymmetric e + e~ stor- 
age ring [8]. This represents an integrated luminosity 
of 425 fb" 1 collected on the T(45) resonance (on-peak), 
which corresponds to (466 ± 5) million BB pairs. In this 
letter, we present an update of our previous measure- 
ments of the branching fraction B and CP asymmetries 
of B° — > J/tpw [9], which had been performed using an 
integrated luminosity of 232 fb _1 . Belle has also studied 
this mode and has published a branching fraction and 
a time-dependent CP violating asymmetry result using 
29.4 fb _1 and 484.3 fb _1 of integrated luminosity, re- 
spectively [10, 11]. 

We reconstruct B° — > J/ip ir° decays from combinations 
of J/ip — > £ + l~ (I = e, fi) and tt° — > 77 candidates. A de- 
tailed description of the charged particle reconstruction 
and identification can be found elsewhere [12]. For the 
J/ip — > e+e~ (^+/i _ ) channel, the invariant mass of the 
lepton pair is required to lie between 3.06 and 3.12 GeV/c 2 
(3.07 and 3.13GeV/c 2 ). Each lepton candidate must be 
consistent with the electron (muon) signature in the de- 
tector. We form it — ► 77 candidates from clusters in 
the electromagnetic calorimeter with an invariant mass, 
to 77 , satisfying 100 < to 77 < 160 MeV/c 2 . These clus- 
ters are required to be isolated from any charged tracks, 
carry a minimum energy of 30 McV, and have a lateral en- 
ergy distribution consistent with that of a photon. Each 
7T° candidate is required to have a minimum energy of 
200 MeV and is constrained to the nominal mass [13]. 



We use two kinematic variables, /ties an d AE, in order 
to isolate the signal: toes = \J ( s /2 + Pi • Pb) 2 /Ef — p 2 B 
is the beam-energy substituted mass and AE = Eg — 
yfs/2 is the difference between the £>-candidate energy 
and the beam energy. Here the B° — > J/ipn° candidate 
(-B rec ) momentum ps and four-momentum of the initial 
state (Ei,pi) are defined in the laboratory frame, E* B is 
the B rcc energy in the center-of-mass (CM) frame, and 
y/s/2 is the beam energy in the CM frame. We require 
toes > 5.2GeV/c 2 and -0.1 < AE < 0.3 GeV. The 
asymmetric AE cut is used in order to reduce background 
from B meson decays to final states including a J/ip me- 
son, where one or more of the particles in the final state 
is not reconstructed as part of B rcc . 

A significant source of background is from e+e~ — ► 
QQ il = u,d,s,c) continuum events. We combine several 
kinematic and topological variables into a Fisher discrim- 
inant {!F) to provide additional separation between signal 
and continuum. The three variables cos(6h), Lq, and L2 
are inputs to J 7 , where Oh is the angle between the pos- 
itively charged lepton and the B candidate momenta in 
the J/ip rest frame. The variables Lq and Li are the 
zcroth- and second-order moments; Lq = X^|p*| and 

L 2 = J2 t |P*I ( 3 cos2 0« - !) A where P* are the CM mo- 
menta of the tracks and neutral calorimeter clusters that 
are not associated with the signal candidate. The 0j are 
the angles between p* and the thrust axis of the sig- 
nal candidate. We use data collected 40 MeV below the 
T(4S I ) resonance to model background from continuum 
events, and signal Monte Carlo (MC) simulated data to 
calculate the coefficients used in T . 

We use multivariate algorithms to identify signatures 
that determine (tag) the flavor of the decay of the other 
B in the event (Stag) to be either a B° or B°. The fla- 
vor tagging algorithm has seven mutually exclusive cate- 
gories of events and is described in detail elsewhere [14]. 
The total effective tagging efficiency of this algorithm is 
given by ^ e;(l - 2ui) 2 = (30.5 ± 0.4)%, where e t is the 
efficiency of a tag, Wj is the probability of mis-identifying 
a tag, and i runs over the seven tag categories. 

The decay rate / + (/_) of neutral decays to a CP eigen- 
state, when £? tag is a B° (B°), is: 

-|At|/r B o 

f±(At) = — [l±S , sin(Am d Ai)=FCcos(Am (i Ai)], 

4t b o 

(1) 

where At is the difference between the proper decay times 
of the B rec and B tag mesons, t b o = 1.530 ± 0.009 ps is 
the B° lifetime and Am d = 0.507 ± 0.005 ps" 1 is the B°- 
B° oscillation angular frequency [13]. The decay width 
difference between the B° mass eigenstates is assumed to 
be zero. 

The time interval At is calculated from the measured 
separation Az between the decay vertices of B rcc and 
£>tag along the collision axis (z). The vertex of £> roc is 
reconstructed from the lepton tracks that come from the 
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J/ip; the vertex of B tag is constructed from tracks in the 
event that do not belong to B Tec , with constraints from 
the beam spot location and the B [cc momentum. We 
accept events with |Af| < 20 ps whose uncertainty a (At) 
is less than 2.5 ps. 

After the selection criteria mentioned above are ap- 
plied, the average number of candidates per event is ap- 
proximately 1.1 in data. The multiple candidates per 
event result from having more than one choice of ir° per 
event, so we choose the one whose value of closest 
to the 7T° mass reported by the PDG [13]. Overall, the 
true signal candidate is correctly identified 99.6% of the 
time for signal MC simulated data. After this step, the 
signal efficiency is 19.3% and a total of 1120 events are 
selected in on-peak data. 

In addition to signal and continuum background 
events, there are also £?i?-associated backgrounds present 
in the data. We consider B backgrounds from the follow- 
ing types of event: (i) B° -> J/i>K% (ii) B° -> J/tPK*°, 
(iii) B ± -> J/ipK^, (iv) B ± -> J/ipp ± , (v) B° -> J/typ , 

(vi) other B decays to final states including a J/tp , and 

(vii) B meson decays to final states including charm 
mesons. The yields of these backgrounds are fixed to ex- 
pectations (16.2, 9.4, 8.8, 2.3, 0.3, 79.4, and 60.4 events, 
respectively), using branching ratios from world aver- 
ages [15]. We allow these to vary in turn when evalu- 
ating systematic uncertainties. Backgrounds from other 
B decays are small, and have been neglected. 

The signal yield, S, and C are simultaneously ex- 
tracted from an unbinned extended maximum-likelihood 
(ML) fit to the on-peak data sample, where the discrim- 
inating variables used in the fit are ttjes, AE, T and 
At. For each candidate-type (signal, continuum, and the 
aforementioned B backgrounds) we construct a probabil- 
ity density function (PDF) that is the product of PDFs 
in each of these variables, assuming that they are un- 
correlated. These combined PDFs are used in the fit to 
the data sample. The continuum-background mES; AE, 
T , and At PDF parameters are floated in the final fit 
to the data. For all other types the PDF parameters are 
extracted from high-statistics MC samples. The mss dis- 
tributions for signal and B° — > J/ip K® events peak at the 
B mass, and are described by a Gaussian with a low side 
exponential tail (GE). The toes PDFs for all other back- 
grounds are described by ARGUS functions [16]. The 
signal AE distribution is described by a sum of a GE 
distribution and a second order polynomial. We use a 
smoothed histogram of MC simulated data to describe 
the AE PDFs for B° -> J/tpK°, B ± -> J/ipp ± , and B 
meson decays to final states including charm mesons, and 
second order polynomials for the AE PDFs of all other 
backgrounds. We parameterize the T distribution for sig- 
nal and continuum events using the sum of a Gaussian 
and a Gaussian with different widths above, and below 
the mean. The T distributions for all other background 
PDFs are Gaussians. The signal At distribution is de- 



scribed by Eq. (1) convolved with three Gaussians (core, 
tail, outliers) which takes into account a (At) from the 
vertex fit, and tagging dilution. The resolution is pa- 
rameterized using a large sample of fully reconstructed 
hadronic B decays [14]. The nominal At distribution for 
the B backgrounds is the same as for signal, except for 
inclusive B and J/i[>K*° backgrounds, which use an ef- 
fective lifetime determined from MC samples of 1.1 ps. 
The continuum background At distribution is described 
by the sum of three Gaussian distributions. The At PDF 
parameters depend on the flavor tag category. The sig- 
nal yield is fitted using known tag efficiencies listed in 
Ref. [14] for each tag category. The continuum yields for 
the seven tagging categories are allowed to vary in the 
ML fit, and the fractions of B background events in each 
category are determined from MC samples. 

After performing tests on the fitting procedure as de- 
scribed in Ref. [17], we fit the data. The results, cor- 
rected for fit bias, are 184 ± 15(stat) signal events, S — 
-1.23±0.21(stat) and C= -0.20 ± 0.19(stat). Figure 1 
shows distributions of jties, AE, and T for the data, 
where the signal is enhanced by selecting AE < 0.1 GeV 
for the mss distribution, and mgs > 5.275 GeV/c 2 for 
the other distributions. These requirements have a rel- 
ative signal efficiency of 98.8% (92.3%) and background 
efficiency of 64% (10.4%) for m ES (AE and J r )_Figure 2 
shows the At distributions for signal B Q and B° tagged 
events. The signal is enhanced by excluding events from 
the tagging category with the largest value of oj, and 
by requiring m E s > 5.275 GeV/c 2 and AE < 0.1 GeV. 
These requirements have a relative efficiency of 70.0% 
(4.4%) for signal (background). The time-dependent de- 
cay rate asymmetry [N(At) -N(At)]/[N(At) +N(At)] 
is also shown, where N (N) is the decay rate for B° (B°) 
tagged events. 




5.2 5.24 5.28 -0.1 0.1 0.2 0.3 -6 -4 -2 2 
m ES (GeV / c 2 ) AE (GeV) Fisher 



FIG. 1: Signal enhanced distribution of (left) mES, (middle) 
AE, and (right) T for the data (points), sum of signal and 
backgrounds (solid line), sum of backgrounds (dashed line), 
and the continuum background (dotted line). 



Table I summarizes the systematic uncertainties on the 
signal yield, S, and C. These include the uncertainty 
due to the PDF parameterization (including the resolu- 
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FIG. 2: The At distribution for a sample of signal enhanced 
events tagged as B° (top) and B° (middle), where dotted lines 
are the sum of backgrounds and solid lines are the sum of 
signal and backgrounds. The time-dependent CP asymmetry 
(see text) is also shown (bottom), where the curve represents 
the measured asymmetry. 



TABLE I: Contributions to the systematic errors on the sig- 
nal yield, S, and C, where the signal yield errors are given 
as number of events. The total systematic uncertainty is the 
quadratic sum of the individual contributions listed. Addi- 
tional systematic uncertainties that are applied only to the 
branching fraction are discussed in the text. 



Contribution 


Yield 


s 


c 


PDF parameterization 


+0.5 
-1.6 


+0.010 
-0.012 


+0.002 
-0.011 


Boost and z-scale 


±1.1 


±0.001 


±0.002 


Beam spot position 




±0.004 


±0.002 


Fit bias 


± 1.5 


±0.021 


±0.014 


B background yields 


±1.2 


±0.029 


±0.013 


CP content of B background 


±0.4 


±0.002 


±0.002 


Tag side interference 




±0.004 


±0.014 


Total 


+2.3 

-2.7 


±0.04 


±0.03 



tion function) , evaluated by varying the signal and back- 
ground PDF parameters within the uncertainties of their 
nominal values. The PDF parameter uncertainties are 
determined from MC samples of signal and background 
events. The uncertainties associated with the Lorentz 
boost, the z-scale of the tracking system, and the event- 
by-event beam spot position are found to be small. We 
determine the fit bias on signal parameters from ensem- 
bles of generated experiments using signal MC simulated 
data, which is generated using the GEANT4-based [18] 
BABAR MC simulation, embedded into MC samples of 
background generated from the likelihood. We apply cor- 
rections to account for the observed fit bias on the signal 
yield, S, and C of -2.7 events, -0.034, and -0.022, re- 
spectively. The uncertainty coming from this correction 



is taken as half of the correction added in quadrature 
with the error on the correction. Most, but not all, of 
the inclusive charmonium final states that dominate the 
inclusive B background are precisely known from previ- 
ous measurements. Their yields are fixed in the fit. As 
a cross check, yields for the B backgrounds are allowed 
to vary one at a time. The sum in quadrature of devi- 
ations from the nominal result is taken as a systematic 
uncertainty. In order to evaluate the uncertainty coming 
from CP violation in the B background, where appro- 
priate, we introduce non-zero S and C for each back- 
ground in turn. The uncertainty due to CP violation in 
B° — > J/ipK® is determined by varying S and C within 
current experimental limits [14, 19]. For B background 
events decaying into final states with charm, we allow 
for a 20% asymmetry, and we allow for 100% asymme- 
tries in all other B backgrounds. We study the possible 
interference between the suppressed b — > ucd amplitude 
with the favored b — » cud amplitude for some tag-side B 
decays [20]. Systematic uncertainties from the effect of 
mis-alignment of the vertex detector and the use of an ef- 
fective lifetime for inclusive B and J/ifjK*° backgrounds 
are found to be negligible. There are additional system- 
atic uncertainties that contribute only to the branching 
fraction. These come from uncertainties for ir° meson re- 
construction efficiency (3%), the J/tp — > branching 
fractions (1.4%), the number of B meson pairs (1.1%), 
and tracking efficiency (1.0%). We apply a correction for 
charged particle identification efficiency (—1.3 ±0.7% for 
J/tp — > e + e~, and — 3.3± 1.0% for J/tp — > decays) 
based on the results of control sample studies using B de- 
cays with J/tjj mesons in the final state. The systematic 
error contribution from MC statistics is negligible. 
We measure 

B = (1.69±0.14(stat)±0.07(syst)) x 10" 5 , 
S = -1.23 ± 0.21 (stat) ± 0.04(syst), 
C = -0.20±0.19(stat)±0.03(syst), 

where the correlation between S and C is 19.7%. We de- 
termine the significance, including systematic uncertain- 
ties, of non-zero values of S and C using ensembles of 
MC simulated experiments as outlined in Ref. [21]. The 
significance of S and C being non-zero is 4. Oct, which 
constitutes evidence for CP violation in B° — > J/i/jtt de- 
cays. The numerical values of S and C are consistent 
with the SM expectations for a tree-dominated b — > ccd 
transition. All results presented here are consistent with 
previous measurements [9-11]. 
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